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A comprehensive study of high-accuracy photoionization cross sections is carried out using the 
relativistic Breit-Pauli R-matrix (BPRM) method for (hi/ + Fe xvn — > Fe XVIII + e) . Owing to its 
importance in high-temperature plasmas the calculations cover a large energy range, particularly the 
myriad photoexciation-of-core (PEC) resonances including the n = 3 levels not heretofore considered. 
The calculations employ a close coupling wave function expansion of 60 levels of the core ion Fe 
XVIII ranging over a wide energy range of nearly 900 eV between the n = 2 and n — 3 levels. Strong 
coupling effects due to dipole transition arrays 2p 5 —¥ 2p 4 (3s, 3d) manifest themselves as large PEC 
resonances throughout this range, and enhance the effective photoionization cross sections orders of 
magnitude above the background. Comparisons with the erstwhile Opacity Project (OP) and other 
previous calculations shows that the currently available cross sections considerably underestimate 
the bound-free cross sections. A level-identification scheme is used for spectroscopic designation of 
the 454 bound fine structure levels of Fe xvn , with n < 10, / < 9, and < J < 8 of even and odd 
parities, obtained using the ab initio BPRM method (compared to 181 LS bound states in the OP 
work). The calculated energies are compared with those available from the National Institute for 
Standards and Technology database, which lists 63 levels with very good agreement. Level-specific 
photoionization cross sections are computed for all levels. In addition, partial cross sections for 
leaving the core ion Fe XVII in the ground state are also obtained. These results should be relevant 
to modeling of astrophysical and laboratory plasma sources requiring (i) photoionization rates, (ii) 
extensive non-local-thermodynamic-equilibrium models, (iii) total unified electron-ion recombination 
rates including radiative and dielectronic recombination, and (iv) plasma opacities. We particularly 
examine PEC and non-PEC resonance strengths and emphasize their expanded role to incorporate 
inner-shell excitations for improved opacities, as shown by the computed monochromatic opacity of 
Fe xvn . 

PACS numbers: 32.80.Fb,96.60. Jw,52.70.-m 



I. INTRODUCTION 

High precision studies of photoionization of atoms in 
a large number of excited levels are of interest in sev- 
eral areas. Derivative quantities such as laboratory and 
astrophysical plasma opacities [lHj], spectral models of 
photoionization dominated soures such as active galactic 
nuclei [f|, non-local-thermodynamic-equilibrium models 
of stellar atmospheres Q , and total electron-ion recom- 
bination rates [7|, depend on the accuracy of underlying 
treatment of photoionization for the entire atom(ion) at 
all energies where it is abundant under specifhc plasma 
conditions. 

There have been a number of theoretical calculations 
to produce large amounts of data for practical appli- 
cations, in particular under the Opacity Project (OP) 
3 and the Iron Project (IP) (32|. But they are usu- 
ally limited in terms of accuracy owing to computational 
constraints, and in detailed examination of the photoab- 
sorption process for excited levels. Furthermore, except 
for light elements or simple few-electron systems, large- 
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scale calculations do not adequately treat the two most 
important atomic effects, relativistic fine structure and 
resonances. Many previous calculations are basically in 
LS coupling 3], or which take account of fine structure 
through an algebraic transformation rather than a rela- 
tivistic calculation per se. The primary reason of course 
is the difficulty of doing so, since the calculations for 
photoionization of a large number of fine structure levels 
remains an enormous task. Close coupling calculations 
that account for the ubiquitous resonances in photoion- 
ization cross sections are mostly done using the R-matrix 
method pj-fiol]. The extension to include relativistic ef- 
fects led to the development of the Breit-Pauli R-matrix 
(BPRM) codes [IJ, [l2| , and to aforementioned work un- 
der the Iron Project for relatively simple atomic systems 
such as the He- like and Li-like ions (e.g. [H[). 

Such large-scale calculations with high precision may 
be of crucial importance in resolving a major astrophys- 
ical problem [14[ that is related to plasma opacities [15| 
and recent studies of solar abundances of common el- 
ements such as C, N, O and Ne. The new abundances 
are widely discordant with the standard solar abundances 
by up to 50% [14]. The new abundances are derived 
from high-resolution spectroscopy and the most advanced 
3D hydrodynamic non-local-thermodynamic-equilibrium 
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(NLTE) models However, such a huge decrease in solar 
elemental composition is not supported by heliosmologi- 
cal observations of solar oscillations measured with great 
accuracy (l6j . In addition, predictions of solar interior 
models are also in conflict with the new abundances [l7| . 
Given the inverse relationship between opacity and abun- 
dances, tt has been suggested that a marginal, but highly 
significant, increase in mean opacities of about 10-20% 
could possible reconcile the new and the standard solar 
abundances [l8j . 

Another need for higher precision opacities arises due 
to the recent capability of creating stellar interior con- 
ditions in the laboratory Inertial confinement fusion 
(ICF) devices can now produce plasmas in LTE at tem- 
peratures and densities approaching those at the bound- 
ary of the radiative zone and the convection zone (CZ) 
in the Sun. Measurements have been carried out at the 
Z-pinch facility at the Sandia National Laboratory of 
monochromatic opacity spectra of iron at T e ~ 150-200 
eV and N e ~ 10 22_23 /cc. Fe xvn is the largest contrib- 
utor to opacity amongst the iron ions prevalent at that 
temperature and density, from Fe xvi - Fe xix with ion- 
ization fractions 0.0.0757, 0.302, 0.366 and 0.18 respec- 
tively. However, there are marked differences between the 
Z-pinch observations and theoretically computed opaci- 
ties [l|. Thus, both astrophysical and laboratory obser- 
vations require high-precision studies to understand or 
resolve the discrepancies. 

We now note one of the basic atomic physics assump- 
tions made in earlier theoretical calculations of plasma 
opacities: resonances are treated as lines, and perturba- 
tive treatments are then employed to include continuum 
contribuions. Depending of the temperature and density, 
most of the opacity could stem from inner-shell excita- 
tions in complex atomic species such as iron ions with 
closed L and M subshells [2, 19]. Such photoexcitation- 
of-core (PEC) process occurs at high energies above the 
first ionization threshold. Therefore, it manfiests itself as 
resonances in the bound-free. The corresponding PEC 
resonances are large features present at incident photon 
frequencies corresponding to dipole transitions in the core 
ion [H, |27|. It follows that the major source of opac- 
ity in any plasma source, astrophysical or in the labora- 
tory, lies in bound-free resonances in general, heretofore 
treated mostly as bound-bound transitions or lines, and 
in PEC resonances in particular. But the approximation 
of resonant photoionization by lines was made in the OP 
0, |i| and OPAL0 projects because opacities calcula- 
tions are extremely intensive, and require huge amounts 
of atomic data for all consituent elements in many exci- 
tation/ionization states. Nevertheless, good agreement is 
found between the two projects for mean opacities over a 
wide range of temperatures and densities ([Hi, discussed 
later). But now, owing to the major problems outlined 
above, it appears necessary to compute atomic parame- 
ters to higher precision to address this issue, albeit on a 
relatively small scale for individiual ions to begin with. 

This report presents pilot and prototypical calculations 



required for revised opacities work. Ne-like Fe XVII is a 
prime contributor to solar opacity near the convection 
zone boundary [l[ , with many PEC features in photoion- 
ization cross sections. By carrying out one of the most 
comprehensive BPRM calculations for photoionization of 
an atomic system, we have computed and analyzed level- 
specific fine structure cross sections, and the distribu- 
tion of resonances over a wide energy range encompassing 
complex atomic structures. In previous studies of pho- 
toionization including resonances, Fe xvn was studied 
mainly in the low energy region using a small wavefunc- 
tion representation. An R-matrix calculation in LS cou- 
pling using a 2-state wave function expansion for the core 
ion Fe xvm was carried out under the Opacity Project 
(M. P. Scott, see [39|). A later coupled channel calcu- 
lation, aimed at unfied electron-ion recombination cal- 
culations for (e + Fe xvm ) — > Fe xvn , employed the 
(BPRM) method using a 3-level (or 3CC) wave function 
expanion by [20j . The recombination spectrum obtained 
using 3CC photoionization compared very well with the 
experimentally observed spectrum in the very low energy 
region. But due to the ion's abundance in high temper- 
ature plasmas, a preliminary calculation using a 60-level 
expansion for Fe xvm was initiated by [2l|. This initial 
study on photoionization cross sections focused only on 
one symmetry: levels with total J=0. But to examine its 
opacity in toto a much more extensive calculation is nec- 
essary. This study is a full-scale photoionization calcula- 
tion using the large 60-level wave function expansion for 
Fe xvm , with energy levels up to approximately 1 keV. 
That energy range should account for nearly all of the 
resonance structures due to L-shell excitation of Fe XVII 
at temperatures ^200 eV where it is abundant near the 
solar CZ; the Planckian blackbody distribution decreases 
to render Fe xvn ionization fraction and opacity to be 
insignifiant beyond E ~ 1.5 keV. We especially aim to elu- 
cidate photoexcitation-of-core (PEC) resonance features 
that greatly affect photoionization in the high energy re- 
gion. 

In this report we point out other potential applications 
as well. Since the cross sections are computed up to levels 
with n(SLJ) < 10, they should be useful for non-LTE 
spectral modeling in the X-ray region. The present work 
yields accurate energies and spectroscopic identification 
of the 454 bound levels of Fe xvn computed through 
ab initio calculations. In addition, level-specific partial 
photoionization cross sections of Fe xvn into the ground 
state of Fe xvm are computed, as required for future 
studies of unified recombination cross sections including 
both radiative and dielectronic recombination processes. 



II. THEORY 

Photoionization calculations have been carried out in 
the close coupling (CC) approximation using the R- 
matrix method as developed under the Opacity Project 
di] and the Iron Project (HH^. In the CC approxima- 
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The first term on the right in Eq. (1) represents both 
bound states and free (continuum) states of the (e + ion) 
system. If all channels are closed then the state is bound 
an represented by ^b', asymptotically, all channel func- 
tions are expnentially decaying. On the other hand, a 
continuum state corresponds to some channels open and 
some closed, and is referred to as typ. Asymptotically, 
the open channel functions are oscillating, as for a free 
electron. The transition matrix element for photoioniza- 
tion is 

< * fl ||D||* F >, (5) 

where the sum is over the ground and excited states where D = J2i r i is tnc dipole operator and the sum is 

of the target or the core ion. The (N+l)" 1 electron over the number of electrons; and ^p are the bound 

with kinetic energy k 2 corresponds to a channel labelled and continuum wave functions. The transition matrix 

SiLi(Ji)Trikf£i(SL(J)Tr). The $ 3 s are bound channel element can be reduced to generalized line strength as 
functions of the (N+l)-clcctron system that account for 
short range correlation not considered in the first term 

and the orthogonality between the continuum and the S = | < ^||Di,||^j > | 2 = 
bound electron orbitals of the target. 

The relativistic Hamiltonian in the Breit-Pauli R- 

matrix (BPRM) approximation is given by where ^ and \Pf are the initial and final state wave func- 
tions. The photoionization cross section (<Jpi) is propor- 

i?jv+i = tional to the generalized line strength as, 

EN+l J V72 1Z , ^N+l 2 1, jrmass _i_ rrDar r rrso 
i=l \ V i ~ 77 "+" l^j>i 7~~ J + U N+l + n N+l "+" 11 N+V 4tt 2 1 

. ( 2 ) apl = 17^ S ' (?) 
where the last three terms are relativistic corrections, 

respectively: where g is the statistical weight factor of the bound state 

and oj is the incident photon energy (Ry). 

the mass correction term, iJ mass = — 2L. J^iPii 
the Darwin term, H Dal = J2i V 2 (f ), 

the spin - orbit interaction term, H°° = Za 2 £ ( ^lj.ss. IIL COMPUTATIONS 

(3) 

Eq.(3) representes the one-body terms of the Breit in- BPRM photoionization cross sections of Fe XVII were 

teraction. Another version of "full" BPRM codes has computed in the 60CC expansion over the core ion Fe 

been developed including the two-body terms, but those xvm > and free-electron wave functions with partial 

are only of importance in forbidden (E2, Ml, M2) tran- waves up to / < 12. The one-electron basis set of orbitals 

sitions where two-electrons correlations may play an im- used to construct a configuration-interaction representa- 

portant role In contrast, plasma opacities m toto tion of the eigenfunctions of Fe xvm were obtained from 

are determined by strong El (dipole allowed and inter- the code SUPERSTRUCTURE 29], which employs a 

combination) considered in this work. Thomas-Fermi-Dirac-Amaldi central field potential. The 

Substitution of * B {e+ion) in the Schrodinger equation choice of the one-electron orbital basis set is not unique; 

Hartree-Fock orbitals of Slater type orbitals are also used 

Hn+i&e = E&p (4) m R-matrix calculations H-E3- At sufficiently high den- 
sities, viz. N e > 10 22 cc, the electronic orbitals may 

introduces a set of coupled equations that are solved us- themselves be altered by plasma effects such as Debye 

ing the R-matrix method. The solution is a continuun screening. However, we do not consider those effects in 

wavefunction "J p for an electron with positive energies (E OP calculations which entail the isolated atom app roxi- 

> 0), or a bound state ^b at a negative total energy (E < mation and the so called "chemical picture" 0, |33| (dis- 

0). The complex resonance structures in photoionization cussed later). 

cross sections result from channel couplings between the Table 1 gives the set of 11 configurations with filled 

continuum channels that are open {kf > 0), and ones Is 2 shell that was optimized to obtain Fe xvin core wave 

that are closed (kf < 0). Resonances occur at electron functions [2l[. It lists the 60 levels included in the 60CC 

energies k 2 corresponding to autoionizing states belong- expansion. For a small number of levels the available en- 

ing to Rydberg series, SiLiixivl where v is the effective ergy level data [26] compiled by the National Institute 

quantum number, converging on to the target threshold for Standards and Technology (NIST), is found to be in 

SiLj. very good agreement with theoretical calculations [2lj . 



tion the atomic system is represented as the 'target' or 
the 'core' ion of N-electrons interacting with the (N+l)" 1 
electron. The (N+l)*' 1 electron may be bound in the 
electron-ion system, or in the electron-ion continuum de- 
pending on its energy to be negative or positive. The 
total wavefunction, *S?e, of the (N+l)-electron system in 
a symmetry Jtt is an expansion over the eigenfunctions 
of the target ion, Xi i n specific state SiLi(Ji)iTi, coupled 
with the (N+l)*' 1 electron function, 6^: 

^ E {e + ion) = A^Xi{ion)6 l + ^c ] <5> :j: (1) 

i 



N+l 



(6) 
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TABLE I: Energy levels of the target ion Fe xvm included 
in the eigenfunction expansion of Fe xvn . Note the large 
energy gap of ~47 Ry between the n = 2 and n = 3 com- 
plexes. The target was optimized with a set of 11 con- 
figurations with open L-shell but closed K-shell: 2s 2 2p J (l), 
2s2p 6 (2), 2s 2 2p 4 3s(3), 2s 2 2p 4 3p(4), 2s 2 2p 4 3d(5), 2s2p 5 3s(6), 
2s2p 5 3p(7), 2s2p 5 3d(8), 2p 6 3s(9), 2p 6 3p(10), 2p 6 3d(ll). 



i Config- Term 2J E(Ry) i Config- Term 2J -E(Ry) 
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One of the main points to note, and which is of consider- 
able significance in the energy distribution of resonances, 
is the clustering of 60 fine structure levels in the n = 2 
and n — 3 complexes, and the large energy gap of ^50 
Ry between them. Although the target energies are quite 
accurate, some further improvement in the positions of 
resonances is also achieved by replacing the Fe XVIII level 
energies with observed ones, wherever available [26| . dur- 
ing diagonalization of the (N-t-l)-electron Hamiltonian of 
the (e + Fe xvm ) system. 

The second sum in the wavefunction expansion given in 
Eq. (1) is the bound-channel term consisting of selected 
electronic configurations for the electron-ion system. We 
include 27 configurations of (A+l) -electron bound chan- 
nels of Fe xvn , specified by a range of mininum and 
maximum occupancies (listed within parentheses after 
the orbitals) as: 2s(0 - 2), 2p(3 - 6), 3s(0 - 2), 3p(0 - 2) 
and 3d(0 — 2). All SLJtt symmetries of the electron-ion 
system formed from the target states coupled with an in- 
teracting electron with continnum partial waves < I < 
12 are considered. 

The ab initio bound state energies of the electron-ion 
system computed by the R-matrix codes in intermediate 



SLJ coupling are not spectroscopically identified a pri- 
ori. For complex ions it is a highly non-trivial task to 
assign LS term and SLJ level designations. It is par- 
ticularly difficult for fine structure levels computed in 
BPRM calculations owing to near-degeneracy of levels in 
high-Z or high-z ions, particularly for high-n and The 
code PRCBPD [37] is employed for level identification, 
using the information computed by the BPRM codes, as 
explained in the next section. In addition, the BPRM 
energies are cross checked against observations and SU- 
PERSTRUCTURE whereever possible. However, we still 
note that this is one of the most laborious tasks requiring 
some judgement in the final assignments for levels highly 
mixed by configuration interaction. 

The present calculations cover a large energy range and 
variations in the distribution of resonance complexes. In 
the near-threshold region, below n — 2 levels, the upi 
were resolved on a fine energy mesh, with 4000 energies 
upto 0.4 Ry above the ionization threshold. However, 
such fine resolution is computationally prohibitive for all 
symmetries and levels over the entire range of ~65 Ry 
where resonances due to L-shell excitation occur. There 
are energy regions where the resonances are sparse, as 
opposed to regions where they are densely clustered. A 
number of energy meshes are used to ensure that the 
overall, as well as the detailed, contribution of resonances 
to opi is taken into account. At photoclcctron energies 
above all 60 target thresholds, the api are extrapolated 
as in [36j . 



IV. RESULTS AND DISCUSSION 

The results of Fe XVII are divided into a few subsections 
below. 



A. Fine Structure Levels and Oscillator Strengths 

of Fe xvn 

Present BPRM calculations for fine structure bound 
levels are intended to form a complete set for most practi- 
cal applications. As explained in an earlier paper |34| . the 
high-lying excited core states do not form bound states of 
the electron-ion system, and corresponding channels have 
insignificant effect on the bound state energies. This is 
particularly true of multiply ionized ions, where there are 
large gaps between the ground complex and the next ex- 
cited complex of levels; the energy separation increases 
as zr. The core ion Fe XVIII has an energy gap of ^47 Ry 
between the n=2 and n=3 levels (Table 1). Although the 
present study includes both the n=2 and 3 complexes, all 
bound levels of Fe xvn have the Fe XVIII parent level 
as the ground level or another n—2 excited level, i.e. 
2sV 2 P? /2<3/2 , 2s2p e 2 S 1/2 . 

R-matrix calculations for bound energy levels entail 
a 'search' for zeroes of the electron-ion Hamiltonian 



5 



32, 38]. where the eigenvalues of the (N+l)-electron sys- 
tem occur as in Eq. (4). A total of 454 eigenenergies of Fc 
XVII were found, subject to the choice of n < 10, < I < 
9, and < J < 8 of even and odd parities. These levels 
have been identified spectroscopically using a numerical 
procedure implemented in the code PRCBPID L 37j. It is 
based on (i) detailed analysis of quantum defects along 
Rydberg series of levels, (ii) parentage of the core ion 
states, (iii) fractional or percentage contribution of closed 
channels which translates into configuration interaction 
of corresponding bound electron-ion configurations, and 
(iv) angular momentum algebra. The level and ener- 
gies are similar to the 3CC case carried [2l|, with some 
differences because of mixed levels and near-degenerate 
quantum defects. 

For relatively few levels the present calculated ener- 
gies of Fe XVII are compared with the observed values 
available in the NIST compilation [26| in Table 2. The 
overall agreement is very good, ^0.1% or better, for most 
levels including ones that are highly excited. The num- 
ber next to the J-value in Tabel 2 (Col. 3) is the rela- 
tive position of the correspoding calculated energy in its 
own symmetry Jit. As an example we discuss levels of 
a particular summetry SLJ i 1 P^. In the present work 
we obtain 35 1 P 1 ° levels, as opposed to 11 given in the 
NIST tables; they are compared in Table 2. Of those 11 
energies, 5 energies agree to ^0.1% or better. That in- 
cludes highly excited levels such as 2s 2 2p 6 4p 1 P°, listed 
as the 19" 1 computed level (Col. 3). For two other levels, 
2s 2 2p 5 6d 1 Pf and 2s 2 2p 5 5d 1 Pf, the energies agree to ~ 
0.2% and 0.3% respectively. However, the energies of the 
remaining three levels differ by 1-4%. The largest dis- 
crepancy is for the level 2s 2 2p 5 5s 1 P 1 °. The complete set 
of 454 energies with spectroscopic designations is avail- 
able electronically (at the NORAD atomic data website 
given at the end of the paper). 

We also report that the oscillator strengths for elec- 
tric dipole transitions among all the bound levels are ob- 
tained and are available electronically. The earlier set 
corresponds to 3-CC calculations [4l| . 



Photoionization Cross Sections of Fe xvn 



TABLE II: Comparison of calculated energies, E c , of Fe XVII 
with the measured values, E [2(|. ij indicates position of the 
calculated level for symmetry J. An asterisk next to a level 
indicates incomplete set of observed levels for the state. 



Conf Term J : i,j 



Sc(Ry) 



2s 2 2p b 
2s 2 2p 5 3s 
2s 2 2p 5 3s 
2s 2 2p 5 3s 
2s 2 2p 5 3s 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p 
2s 2 2p 5 3p X D 
2s 2 2p 5 3p 1 S 



3po 

3 r 

3po 
Ipo 

3 S 
3 D 
3 D 
3 D 
!p 

3p 
3p 
3p 
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Photoionization cross sections of all 454 bound fine 
structure levels of Fe xvn are computed and analyzed 
in a variety of ways, in particular the distribution of res- 
onances that lie in the large energy gap between the n 
= 2 and n = 3 complexes. This is a comprehensive set 
of cross sections computed using the BPRM method, as 
required for applications such as plasma opacities, syn- 
thetic spectral models, and spectral diagnostics of X-ray 
absorption lines observed from astronomical objects. The 
BPRM calculations yield about two and a half times as 
many levels as the previous OP calculations in LS cou- 
pling [3{|, which resulted in 181 LS bound states with 
n(SL) < 10. Some important charateristic features of 
the BPRM cross sections are illustrated and discussed 



below. 



1. Photoionization of the Ground State 

First we consider the ground state of Fe xvn . Fig. 1 
presents photoionization cross sections opi of the ground 
state 2s 2 2p 6 1 S. While panel (a) presents the cross sec- 
tion in LS coupling under the Opacity Project (M. P. 
Scott, see TOPbase data [3{|), panels (b,c) present apj 
from the relativistic BPRM calculations; (b) presents the 
total cross section, and (c) presents partial cross section 
for leaving the core ion Fe xvni in the ground state 
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FIG. 1: Photoionization cross sections, api, of the ground 
level ls 2 2s 2 2p 6 ( L So) of Fe xvn : (a) total a PI in LS cou- 
pling from the Opacity Project database TOPbase [391 ] - (b) 
present total opi in the relativistic BPRM formulation, (c) 
present partial api with residual ion Fe xvin in the ground 
state. Present calculations including relativistic fine structure 
provide more accurate energies and more highly resolved res- 
onances. The background for the partial cross sections (c) is 
lower that the total (b) because the channels for leaving the 
core ion in various excited levels are excluded. 



following photoionization. The earlier total opi from 
the two-state R-matrix OP calculation, and the present 
BPRM using a 60-level expansion are of similar magni- 
tude for the background. But due to a much larger num- 
ber of fine structure channels in the BPRM calculation, 
panel (b) shows many more resonances within the n=2 
complex (as also found in the earlier 3CC calculation). 
Also, a much finer energy mesh in the present work has 
resolved the resonances more completely. Resolution of 
resonances in the low energy region, especially at and 
near threshold is crucial for calculating recombination 
and photoionization rates. The ground state opi at high 
energies decreases slowly, showing an insignificant effect 
due to highly excited n=3 core states in the 60CC calcu- 
lation, except for introducing small and weak resonance 
structures. 

The partial cross sections in panel (c) corresponds to 



photoionization of Fe xvn leaving the ionized core Fe 
xvin only in the ground state. The background is lower 
than that of the total apj in (b) because there are no 
additional contributions from photoionization channels 
into excited levels. However the resonance features are 
the same for the total and the partial cross sections until 
the first excited state 2s2p 6 1 S Q of Fe xvin . 

2. Photoionization of Excited Levels 

Although the background ground state cross sections 
remain about the same between the 3CC and the 60CC 
calculations, the excited states are considerably more af- 
fected. One of the main results of the present BPRM 
calculations with a large wavefunction expansion is that 
photon absorptions and core excitations to the n=3 lev- 
els produce extensive resonance structures, in contrast to 
earlier calculations. This contradicts the assumption that 
channels of highly excited core states, especially when an 
energy gap like in the present case exists, might be too 
weakly interacting to produce any significant effect. 

Fig. 2 presents total photoionization cross sections of 
excited levels 2s 2 2p 5 3p( 1 P 1 °) (a,b), and 2s 2 2p 5 3d( 1 D%) 
(c,d). These two levels are chosen because they have a 
single fine structure component, so their cross sections 
can be directly compared with the earlier OP results for 
the LS terms 2s 2 2p 5 3p( 1 P°) and 2s 2 2p 5 3d( 1 D°). The 
top panels (a,c) of Fig. 2 present the OP results, and the 
bottom panels (b,d) from the present work. Both of these 
excited levels show resonances due to core excitations to 
the n — 2 and n = 3 levels that are not present in the 
OP results. The features reveal how crucially important 
the set of resoances belonging to n=3 core thresholds are, 
compared to those of n = 2. The resonance peaks are or- 
ders of magnitude higher than those from the n=2 series. 
There is only a relatively small gap with smooth back- 
ground without resonances, between the highest thresh- 
old of the n=2 levels and the appearance of lowest of the 
n=3 complexes of resonances. While the background is 
enhanced considerably in the BPRM results by strong 
resonances, the featureless and slowly decreasing cross 
sections from the OP miss out all those features and 
vastly underestimate opi. 

Resonances in photoionization cross sections corre- 
spond to Rydberg series of autoionizing levels at energies 
E res — E c vl, where E c is an excited core threshold and Vi 
are the efeective quantum number and angulur momen- 
tun of the interacting electron. The Rydberg resonances 
are narrow and lie below the threshold E Cl approximately 
at energies given by the simple expression 

E res = E c -z 2 /v 2 . (8) 

However, strong PEC resonances manifest themselves in 
apj of excited bound levels of the electron-ion system 
[271 [35j . The PEC resonances are wide and occur at 
energies where the ground state of the core ion undergoes 
a strong dipole allowed transition. Among the 60 target 
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FIG. 2: Comparison of present photoionization cross sections 
a pi with those from the OP (see [39() for excited singlet lev- 
els (no fine structure): (a,b) 2p 5 3p 1 P, and (c,d) 2p 5 3d( 1 D°). 
In contrast to api from OP, the present results demonstrate 
that without inclusion of n=3 core levels of Fe xviii the cross 
sections are considerably underestimated throughout most of 
the energy region of interest in practical applications. Reso- 
nances are included as lines in the OP work [3J, as in other 
opacities calculations [!]• 



or core levels for Fe xviii there are 30 such transitions as 
given in Table 3. Hence, in photoionization cross sections 
of Fe xvii there are 30 possible PEC resonances, most 
with overlapping profiles. 

Fig. 3 presents photoionization cross sections of three 
highly excited levels of Fe xvii and demonstrably large 
PEC features, 2s 2 2p 5 nf( 3 D 1 ) with nf = 5f, 7f, 9f. Be- 
yond the resonances due to n — 2 core levels, the back- 
ground decreases smoothly for the three levels. However, 
as the resonances due to n = 3 core levels appear, the 
backgound rises and the PEC resonances manifest them- 
selves as prominently high and wide structures. These 
enhancements are related to the radiative decay rates of 
the dipole allowed levels. Table 3 shows most of the de- 
cay rates from n = 3 levels are much larger than those 
from the n = 2 levels, by one or two orders of magni- 
tude, and result in stronger resonances. With hundreds 
of bound levels computed in the present work, the cross 



TABLE III: Dipole allowed and intercombination El transi- 
tions from the ground level 2s 2 2p 5 ( 2 P^ 2 ) to excited states of 
the core ion Fe xviii , and corresponding oscillator strengths 
(/). These transitions introduce the PEC resonances in pho- 
toionization cross sections of Fe xvii . A number 'a(n)' means 
ax 10". The level indices correspond to those of Table 1. 
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sections in the high energy region from about 57 Ry to 
65 Ry is greatly enhanced, mainly by complexes of PEC 
resonances corresponding to core excitations via dipole 
transitions listed in Table 3. 

Of particular importance is the characterstic shape of 
PEC resonances. The resonances encompass an energy 
range of ~10 Ry or well over 100 eV. The shape is de- 
termined by channel coupling effects. In the case of the 
simple limit of an isolated resonance it is like the typ- 
ical Fano profile. However, the PECs are generally af- 
fected by strong coupling among many target levels and 
channels, with a large number of superimposed non-PEC 
resonances converging on to the target levels. The dis- 
tribution of the continuum or the differential oscillator 
strength would not be generally reproduced by an iso- 
lated resonance approximation, and requires a coupled 
channel calculation. 



C. Resonance oscillator strengths 

A fundamental approximation made in existing opac- 
ities calculations is to treat autoionizing resonances as 
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FIG. 3: Total photoionization cross section opi of highly ex- 
cited levels of Fe xvn : 2s 2 2p 5 n/( 3 L>i) with nf = 5/, If, 9/. 
The cross sections decrease monotonically at lower energies, 
but rise considerably while approaching the PEC resonances. 
The figure shows that the cross sections would be enhanced 
all throughout the energy range, compared to currently avail- 
able near-hydrogenic cross sections that decrease much faster 
as an ~ E~ 3 . The largest PEC resonances are marked by 
arrows. Although only two such PECs are pointed out, the 
resonance structures and the resulting enhancement represent 
the combined effect of all PEC resonances in Table 3. It may 
be also be noted that the PEC positions remain the same for 
all Rydberg levels. 



lines. In other words, inner-shell excitations leading 
to bound-free autoionizing states are treated as bound- 
bound transitions. The final state may be further cou- 
pled to a feature-less continuum perturbatively to ob- 
tain autoionization widths at a single energy associated 
with the bound-bound transition. This indepndent res- 
onance approximation neglects the intricate coupling ef- 
fects that are otherwise included via the close coupling 
method. Since most of the contribution to opacities orig- 
inates from inner-shell excitations, with final levels as au- 
toionizing states, their impact on opacities bears closer 
inspection and is briefly discussed below. 

In Fig. 4 we examine in some detail one of the cross sec- 
tions shown in Fig. 3, corresponding to 2s 2 2p 5 7/( 3 £>i). 



FIG. 4: L-shell PEC resonances in photoionization of Fe xvn 
. The resonances cover all 60 ionization thresholds included in 
the CC calculation, including all PEC resonances correspond- 
ing to EI transitions (dipole allowed and intercombination) in 
Table 3. 

The cross sections are delineated at 24,653 energies; the 
energy intervals are chosen so as to resolve resonance 
profiles in so far as practical. Postions of all of the PEC 
resonances due to transitions in Table 3 are marked in 
Fig. 4 as arrows. In addition, there are a large number of 
non-PEC Rydberg resonances converging on to the other 
excited levels. Whereas the few An(2 — 2) resonances in 
Fig. 4a are inconsequential, the myriad An(3 — 2) PECs 
in Fig. 4b and 4c dominate the distribution of oscillator 
strength over a large region ^56-66 Ryd, or >100 eV, due 
to the L-shell excitation array 2p — > 3d. The cumulative 
resonance (or bound-free continuum) oscillator strength 
corresponding to F ig. 4 is related to the photoionization 
cross section (viz. [22]) as 

MCD, — > eSLJ : 3 (P, D, F) (J = 0, 1, 2)°] = (9) 

— i-^ / ° a PI (2s 2 2p 5 7f ^D^de, 
±-Kaa 2 \ Jo 

where e is the energy relative to ionization threshold and 
up to e = ~80 Ry. An integration over the range shown 
in Fig. 4 yields the partial resonance oscillator strength 
including all of the PECs and the non-PEC resonances 
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due to coupling of all 60 Fe xvm levels. The sum over the 
oscillator strengths corresponding to the 30 PECs gives 
a total fpEC = 2.31. The integrated resonance oscillator 
strength f r is found to be 4.38. In other words, the 30 
PEC resonances involving transitions up to the n = 3 
levels of the core ion Fe xvm contribute over half of all 
the continuum bound-free oscillator strength in photoion- 
ization of any excited state of Fe xvn . We note that, 
without loss of generality, we chose an excited level to 
demonstrate the quantitative effect of PEC resonances; 
they manifest themselves in photoionization of most lev- 
els. 



D. Monochromatic Opacities 

Opacity calculations are a complex undertaking that 
requires atomic data for a large number of ions in varying 
plasma environments 0, H[ • One of the main reasons for 
the large-scale calculations carried out under the present 
study is to enable test calculations to benchmark avail- 
able opacities for Fe ions. To that end, the present BPRM 
calculations employ similar cut-offs, n < 10 and £ < n—1, 
as the OP work [3j. However, one usually expects only 
the ground state and low-lying metastable states to be 
significantly populated. The ionization fractions and 
level populations are computed using an equation-of- 
state, such as the modified Boltzmann-Saha formulation 
in the 'chemical picture' [33j], based on the premise that 
isolated atoms exist, albeit perturbed by the plasma en- 
vironment • At low densities and temperatures the ion 
fractions and occupation probabilties of high-lying levels 
are several orders of magnitude smaller than those for 
the ground state and metastable levels. But in the high 
temperature-density regime N e > 10 24 cm" 3 , T e > 10 6 
K, approaching those in stellar cores, electron-ion recom- 
bination rates can be large, and increase rapidly as N 2 . 
Even a small population in excited levels would then be 
susceptible to the resonant enhancements due to PEC 
resonances, which are currently neglected and the cross 
sections for excited levels are taken to be nearly hydro- 
genic, instead of the accurate form exemplified in Figs. 2 
and 3. 

While we have discussed integrated resonances oscil- 
lator strengths embedded in bound-free cross sections, 
there is no direct equivalence or 1-1 correspondence with 
bound-bound oscillator strengths, as generally computed 
in opacities calculations. Among the factors that distin- 
guish the two are overlapping profiles and large energy 
widths of PEC resonances, reflecting the coupling of con- 
tinua belonging to many target levels and strong dipole 
moments among them. Owing to the huge scale of data 
needed, definitive checks can only be made by calculating 
revised opacities using atomic data, as outlined herein. 

A practical problem likely to be encountered is the 
high energy resolution needed to represent resonance pro- 
files. Whereas scattering cross sections are bounded by 
the unitarity condition (viz. [23| . no such upper bound 



exists for individual values of photoionization cross sec- 
tions. The peak photoionization values may rise arbi- 
trarily high, and numerical integration would tend to be 
inaccurate as the resonance profile approaches a Delta 
function; the integral is finite but the width is extremely 
narrow and impractical to resolve. However, the statisti- 
cal methodology adopted in opacities work is to employ 
the opacity sampling technique. Monochromatic opacity 
spectra are sampled at approximately 10,000 points (viz. 
3]), although the atomic data are much more finely re- 
solved. It has been verified that the statistical averages of 
the most important quantity, the Rosseland Mean Opac- 
ity (RMO), do not deviate by more than 1-3% even if the 
atomic cross sections are 'sampled' at 10 5 or 10 6 points 
0. Therefore, the energy mesh of ~30,000 points used in 
this work, predominantly in the region of covered by the 
high energy n — 3 resonances, should suffice for accurate 
opacities calculations. 

Fig 5 presents the monochromatic opacity K(Fe xvn) 
computed using all of the bound-free data in the present 
work, 454 photoionization cross sections resolved as men- 
tioned above. The calculations are carried out us- 
ing a newly developed code for high-precision opacities, 
adapted from the earlier OP work @ (the OP opacities 
are available from the on-line server called OPserver at 
the Ohio Supercomputer Center [11]). The new code 
also employs a frequency mesh of 10 5 points, an order 
of magnitude finer mesh than OP or OPAL, thus obviat- 
ing some resolution issues in the monochromatic opacity 
spectra [39(. Essential components of opacities calcula- 
tions related to bound-bound transitions are retained as 
in the OP code: electron-impact and Stark broadening, 
free-free scattering, and electron-photon interaction in 
the Rayleigh, Thomson, and Compton scattering limits. 
However, resonance profiles have not yet been broadened. 
Whereas bound-bound line shapes are symmetrical, the 
giant PEC resonance profiles span hundreds of eVs (viz. 
Fig. 3 and 4) and are asymmetric. Theoretical formula- 
tions for broadening of resonances and algorithms are still 
being developed. However, as mentioned, plasma effects 
in the chemical picture [33j are included a posteori, and 
do not affect the accuray of the atomic parameters com- 
puted herein. The details of the new opacities code and 
results will be reported later. But we note in passing that 
since the resonances would dissolve more readily than 
lines, it is likely to significantly affect continuum lower- 
ing, manfiest in opacity spectra, and highly senstiive to 
temperature and density. Lines dissolve and eventually 
merge into the continuum in high temperature-density 
regime. Then a precise accounting may not be neces- 
sary; they are often treated as " unresolved transition ar- 
rays" (UTA) since J- J' transitions are merged together 
and often subsumed by line broadening [42j. However, 
our aim is to focus on delineation of atomic features as 
fully as possible so that their contributions to opacity in 
different tcmeperature-density regimes can be accurately 
ascertained. 

The BPRM opacity cross sections in Fig. 5 are com- 
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BPRM With Resonances 




1 i— OP Without Resonances 





Photon Energy E(keV) 



FIG. 5: Partial monochromatic opacity of Fe xvn : logiQK, 
(Mb) at temperature T = 2.24 xlO 6 K and electron density 
7V e = 10 23 cm~ 3 , corresponding to the base of the solar con- 
vection zone where the Fe xvn ion is the largest contributor to 
opacity. The top panel presents results using BPRM bound- 
free cross sections with resonances computed in this work, 
and the middle panel without resonances using earlier data 
from the Opacity Project (OP). Identical datasets from OP 
are employed for the bound-bound transitions, therefore the 
differences are mainly due to the resonances included in the 
present work. For comparison, the bottom panel shows total 
OP opacities [39l ] that include core-excitation resonances as 
lines, with autoionization widths considered perturbatively, 
as well as the high-energy K-shell continuum opacity not yet 
included in BPRM calculations (see text). 



pared to two other sets of results from OP. Although 
over 20,000 oscillator strengths for bound-bound transi- 
tions are also computed, we use exactly the same data 
for lines as the earlier OP calculations described in 
Thus the differences between the the 60CC BPRM results 
from the present work in the upper panel of Fig. 5, and 
the limited OP without resonances in the middle panel, 
are entirely due to the differences in the bound-free data 
sets. The opacity calculations are done at a temperature- 
density representative of the plasma conditions at the 
base of the solar convection zone: logio T (K) = 6.35 
and logio N e = 23. These parameters also lie in the range 
currently under investigation in the Z-pinch plasmas for 



measurment of transmission spectra [l|. In addition to 
resonance contributions, there are some other differences. 
The OP data include extrapolated cross sections out to 
very high energies, ~ 500 Ry. The BPRM data have also 
been processed to include these high-energy "tails", but 
the form is slightly different. The background opacity is 
important to obtain a value for RMOs that spans over 4 
decades in temperature, logio T(K) = 3.5 - 7.5. 

Calculations for the monochromatic opacity k v are car- 
ried out for each ion along isotherms in logio T for a 
range of electron densities logio N e . The Rosseland Mean 
Opacity (RMO) kr is defined in terms of Fly SiS 



1 

KR 



Jo* ' a{u)i: du 

g(u)du 



u 4 e- u 



(1-e-")- 2 , (10) 



where g(u) is the Planck weighting function (corrected for 
stimulated emission). The k v is primarily a function of 
the oscillator strengths /, photoionization cross sections 
<t„, level populations 2Vj, and the line profile factor <j> v , 



—) Nifijfa, K b J = N i( j v . (11) 
m e cj 



Whereas further code developments are required to in- 
clude all opacity contributions, we already find large en- 
hancements in kr due to resonances, primarily from the 
n = 3 complex. The resonance contribution is included 
in existing opacifies codes as inner-shell bound-bound 
transitions [40]. The BPRM value of kr, also including 
the bound-bound oscillator strengths computed in this 
work (as opposed to OP) yields a value of 223.8 cm 2 /g. 
Using the same bound-bound data as OP, the BPRM k- 
value is still 200.3 cm 2 /g, compared to the OP value of 
109.7 cm 2 /g (Fig. 5). The bottom panel in Fig. 5 is the 
total monochromatic opacity spectrum of Fe xvn from 
OP, including all contributions, with the total RMO kr 
= 306.9 cm 2 /g. Thus the BPRM value using the data 
computed in this work is 27% lower. This is primarily 
because of two factors: (I) The high-energy "tails" are 
to made more precise with even more extended BPRM 
calculations including the K-ionization thresholds, which 
would attenuate the bound-free opacity in the large en- 
ergy range between the L and the K levels. Although the 
K-shell opacity is not too significant, the K-shell reso- 
nances would thereby be included. (II) The bound-bound 
oscillator strengths are computed up to n = 10. There- 
fore, there is a relatively small gap in cross sections in 
the region between n = 10 - oo, below the photoion- 
ization thresholds of the 454 bound levels. In the earlier 
OP work, we began the tabulation of cross sections at E 
= —z 2 /v 2 (y ~ 10) below each threshold. Employing a 
similar approximation increases the BPRM RMO value 
from 223.8 cm 2 /g to 260.7 cm 2 /g, to within 20% of the 
total OP RMO. Owing to the signficance of (I) and (II), 
we plan to carry out more extensive calculations than 
in earlier works, as well investigate if higher-n resonance 
complexes with n = 4 or 5 might be needed. 
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E. Unified Electron-Ion Recombination 

One of the important charateristics of PEC resonances 
is that they entail excited states with valence electrons 
that are weakly bound to the core ion. Therefore, dur- 
ing core excitation the outer electron remains essentially 
a 'spectator', temporarily attached to the excited core 
level, and autoionizes as that level decays to the ground 
state. The analogy and connection between the PEC 
resonances and the dielectronic recombination process is 
well known [23|, [53, It is expected that the total 

electron-ion recombination rates of Fe xvn will also be 
commcnsurately enhanced by inclusion of the PEC reso- 
nances via detailed balance (Milne relation, e.g. (23l.[36|). 

A major application of the computed cross sections is 
in benchmarking total unified recombination rate coef- 
ficients, including radiative and dielectronic recombina- 
tion, for Fe xvn at high temperatures, including those 
where Fe XVII is abundant in coronal plasmas, such as 
in the solar corona and solar flares, up to at least 10 7 K. 
That is possible because we have considered a large en- 
ergy range up to the excitation of n = 3 levels of the 
recombining ion Fe xvni . Towards that end, we have 
also repeated the entire photoionization calculation to 
obtain partial cross sections for the 454 bound levels of 
Fe xvn with the residual ion in the ground state alone. 



V. CONCLUSION 

In this report we have presented results from a pilot 
project of complete BPRM calculations for photoioniza- 
tion of an atomic system, Fe XVII, larger than a He- or 
Li-like ion, (hi/ + Fe XVII -> e + Fe XVIII). The aim 
was to study in detail the extent and range of high en- 
ergy resonances of importance in practical applications. 
The BPRM calculations consider all fine structure lev- 
els up to n(SLJ) < 10, with spectroscopic identifica- 
tion. In addition to photoionization cross sections, the 
bound-bound oscillator strengths for transitions among 
the 454 computed levels of Fe xvn are also being com- 
puted. These datasets of atomic parameters should be of 
unprecedented accuracy and generally useful. 

The comprehensive calculations were carried out so 
as to compare with the erstwhile OP calculations that 
treated resonances as lines. A clear distinction is made 
between pure bound-bound transitions and bound-free 
transitions into autoionizing levels. The results from this 
work demonstrate that opacities may be computed using 
BPRM cross sections and transition probabilities, and it 
is likely that plasma opacities in general, and those of 
Fe ions in particular, would be different from earlier ones 
using atomic cross sections that accurately consider the 
energy distribution of resonances. 

Mnochromatic opacities of Fe XVII are computed and 
compared with the OP work. These are not the final 
opaciites, and some further developments are still neces- 
sary, such as resonance broadening mechanisms and K- 



shell contributions. Nevertheless, we are able to compute 
complete BPRM radiative datasets and Rosseland mean 
opacities that are sufficiently close to OP values to im- 
ply that future calculations of plasma opaciites can be 
carried out with higher precision. 



Three other applications of the present work might be 
pointed out. (i) Calculation of accurate unified electron- 
ion recombination rate coefficients, (ii) Benchmarking 
experimental measurements of absolute photoionization 
cross sections on accelerator based light sources, that are 
now being made for multiply charged Fe ions (e.g. 



31 



It has now been established that the experimental beams 
contain ions not only in the ground state but also in sev- 
eral mctastable levels. An admixture of ground state plus 
a few excited metastable levels is therefore necessary to 
benchmark experimental measurements against theory. 
In addition, the experimental measurements are capable 
of reaching energies where PEC resonances occur (e.g. 
[30l| ) . (iii) The radiative data for Fe xvn should be use- 
ful for X-ray spectral diagnostics and non-LTE models 
((Ref. [28[ presents a detailed Grotrian diagram of Fe 
XVII levels up to n=4). 



Estimates of uncertainties in the large sets of param- 
eters reported in this paper, as well as that required in 
future calculations of opacities, are as follows. The ac- 
curacy of the 454 theoretically computed energy levels 
has been ascertained in the text by comparison with ex- 
perimentally observed levels, in general better than 1%. 
The differences between BPRM cross sections and mea- 
surements for photoionization and recombination for Fe 
XVII are found to be within experimental uncertainties 
of ~10-20% [2lJ. Similarly, most of the strong transition 
in the extensive dataset of Fe xvn oscillator strengths 
should be better than 10% accuracy 4l|. The present 
calculations for photoionization employ a much larger 
60CC eigenfunction expansion than the 3CC calcula- 
tions in previous studies [U, Hl[ , and are expected to be 
more accurate. Uncertainites in the BPRM monochro- 
matic opacites should be commensurate with those in 
the underlying 60CC data. However, the final accuracy 
of derived mean Rosseland and Planck opacities ions is 
as yet undetermined since photoionization cross sections 
extrapolated or computed at high energies are to be in- 
cluded for all 454 bound states. The problem is further 
complicated since the aim is to obtain statistical averages 
with sampled opacities (albeit on a finer mesh of 10 5 fre- 
quencies as opposed to 10 in earlier calculations) that 
would approach the accuracy, needed to improve 

over existing uncertainties of ^5% [19j . 



Electronic files for photoionization cross sections, en- 
ergy levels, and oscillator strengths are available electron- 
ically from the NORAD website: www.astronomy.ohio- 
state.edu/^nahar/nahar_radiativeatomicdata 
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